Shear thickening in non-Brownian suspensions: an excluded volume effect 
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Shear-thickening appears as an increase of the viscosity of a dense suspension with the shear rate, 
sometimes sudden and violent at high volume fraction. Its origin for non-colloidal suspension with 
non negligible inertial effects is still debated. Here we consider a simple shear flow and demonstrate 
that fluid inertia causes a strong microstructure anisotropy that results in the formation of a wake 
region with no relative flux of particles. We show that shear-thickening at finite inertia can be 
explained as an increase of the effective volume fraction when considering the dynamically excluded 
volume due to these wake regions. 



The field of complex fluids is diverse and rapidly devel- 
oping with potential for numerous relevant applications. 
Among complex fluids, on one hand we have colloidal sus- 
pensions where Brownian effects play an irnportant role 
while inertial effects are negligible, see e.g. On the 

other hand we have suspensions made out of larger parti- 
cles, (particle radius a > 10/im) where Brownian effects 
are negligible while inertia plays an important role. To 
be specific we shall call this second class of suspensions as 
non-Brownian suspensions or inertial suspensions. Their 
rheology is the topic of this letter. 

Understanding the rheological properties of non- 
Brownian suspensions is not only a challenge from a the- 
oretical point of view Q but has also a significant im- 
pact in many industrial applications, e.g. oil processing, 
cement or coal slurries 0,l8[. 

In one of the earliest work in this field, Einstein showed 
that for a dilute suspension of rigid particles in a New- 
tonian fluid with negligible inertia the relative increase 
in effective viscosity is ^ {5/2)(j) where (p is the volume 
fraction occupied by the particles [see e.g., 0, chapter 
4.11]. For higher concentrations the problem is still not 
well understood. Non-Brownian suspensions may show 
shear-thickenin g, i .e. increase of effective viscosity with 



the shear-rate [lO|, [ll|. If the volume fraction is high 
enough, yet below the geometrical maximum packing 
fraction, 0„j = 0.58 — 0.63, the increase of viscosity with 
shear-rate can be abrupt 1^, i.e. discontinuous shear- 
thickening. 

In this letter, we report on three dimensional Direct 
Numerical Simulations (DNS) of a plane-Couette flow of 
neutrally-buoyant rigid spheres in a fluid. The rheol- 
ogy is governed by two parameters: the volume fraction, 

we use a 



Hi 



(/>, and the shear rate 7. Following Ref. 
non-dimensional form of the shear rate given by the par- 
ticles Reynolds number. Re = pja'^/^xa, where /io, P are 
the fluid viscosity and density and a is the particle ra- 
dius. The effective viscosity is thus function of and 
Re, fj. — fj.Q f((j), Re). For the conflgurations investigated 
here, the effective viscosity, reported in Fig [1] increases 



as the relative strength of the inertial effects (measured 
by Re) increases; a phenomenon we call inertial shear- 
thickening. 

Flow arounds a sphere has been studied exten- 
sively [see e.g., [l3|: for Re ^ 1, Stokes flow, the flow 
velocity is steady and obeys all the symmetries of the 
problem; for Re ~ 1 the flow remains steady but a wake 
forms behind the sphere; as Re is increased further the 
flow becomes unsteady and eventually chaotic. Hence we 
expect that in non-Brownian suspensions with finite in- 
ertia the asymmetry induced by the wake will affect the 
suspension rheology. Indeed, we show that in the wake 
of a first particle there is a vanishing probability to find 
another particle. We obtain an estimate of the average 
volume of the wake region in the suspension by calculat- 
ing the pair-distribution- function. Fig [21 and the relative 
fiux of a pair of spheres. Fig [31 We interpret the vol- 
ume occupied by the wake as an increase in the effective 
volume fraction; this allows us to collapse the data for 
p/fj-o for two different values of (j) into one single func- 
tion of the effective volume (j)c, FigHl This function is 
well approximated by the well-known Filers fit [ll| that 
is an empirical formula describing the variation of viscos- 
ity of a suspension with the volume fraction for vanishing 
inertia. 
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with B = 1.25 - 1.5 and ^ = 0.58 - 0.63 the maxi- 
mum packing fraction l3"la|- Similar data-collapse in 
granular systems under different conditions have been re- 
cently obtained from experimental jl7] and numerical 18 1 
data. We go beyond these studies addressing the prob- 
lem from a microscopical point of view and showing that 
this increase of the effective volume fraction is due to for- 
mation of anisotropic micro-structures characterized by 
the angle-dependent pair-distribution function and mean 
relative particle flux. Note that the existence of such 
micro-structures cannot be inferred from isotropic, angle- 
averaged, observables. Recent investigations [Sj-lsj, |l9| 
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FIG. 1: Normalized effective viscosity /x//io vs particle 
Reynolds number Re\ red solid line with □ symbols </) = 0.21, 
dotted blue line with x symbols cj> — 0.315. 

have stressed the important role played by "hydroclus- 
ters" in shear-thickening in Brownian (colloidal) suspen- 
sions. Here we elucidate the role of the particle clusters 
and microstructure in shear-thickening of non-Brownian 
suspensions with finite inertia. 

We numerically simulate a suspension of rigid spheres 
suspended in a fluid phase described by the incompress- 
ible Navier-Stokes equation. These are solved on a Carte- 
sian mesh in a rectangular box of size 16 a x 16 a x 10 a 
along the spanwise, streamwise, and vertical directions 
{x,y,z), with 8 grid points per a. The fluid is sheared 
in the y — z plane by imposing a constant streamwise 
velocity of opposite sign {l/1)Uo = jh, [2h — 10 a) at 
the two horizontal walls (z = ±/i). Periodic boundary 
conditions are imposed on the other two directions. A 
Lagrangian algorithm is used to solve for the linear and 
angular momentum of the spheres. We impose no slip 
boundary condition on the fluid at the particle surface 
using the Immersed Boundary Method (IBM). Lubrica- 
tion and collision models are employed to capture the 
interaction between spheres when the distance between 
the surface of neighboring particles become smaller than 
the mesh size. The surface of each sphere is discretized 
by about 800 Lagrangian grid points. The code was fully 
validated against several classic test cases, see [2l| for 
more details. Two different values of the volume fraction 
(j) — 0.21 and = 0.315, and several particle Reynolds 
numbers in the range 0.1 to 10 are simulated. 

In Fig. [1] we plot the effective viscosity of the suspen- 
sion, /I, measured as the ratio between the tangential 
stress at the walls and the shear rate 7, as a function of 
Re for the two different values of (p. The effective viscos- 
ity increases with the shear rate (shear-thickens) at fixed 
volume fraction; also it increases with the volume frac- 
tion (j) at fixed Re. Our results are consistent with data 
from recent experimental (20l | and numerical [3] works. 

Next we show that shear-thickening can be interpreted 
as an excluded volume effect. We first calculate the pair 




FIG. 2: Projection of the normalized angle-dependent pair- 
distribution function, g[r ~ 2a, f), on the wall-parallel plane 
(with mean flow from right to left) x — y plane for Re = 0.1 
(a), and Re = 10 (b) for = 0.315. 

distribution function g(r^ f ) that is the probability to find 
a particle pair at given distance r and direction f normal- 
ized by the value for a random arrangement 0, 0, [l^l • 

In Fig. [5] we display g{r,Y) at contact, r — 2a, in the 
wall-parallel x ~ y plane (relative motion from right to 
left), for two different values of Re at ^ = 0.315 (similar 
behavior is observed at lower concentrations). The con- 
tours show that g{r, f ) is not isotropic and the anistropy 
increases as inertial effects, measured by Re, become 
more important. In particular there exists a small wake 
behind the particle where there is a lower probability to 
find a second particle. Increasing the Reynolds number, 
the wake region becomes larger. This depleted region 
acts as an increase in effective volume fraction: this is 
the geometrical volume occupied by the particles plus 
the volume of the wakes (the wake is actually a property 




FIG. 3: Contour plot of particle pair relative flux | q |, Eq. ([2]) 
in the plane of shear y — z for (j> = 0.315: Re = 0.1 (a). 
Re = 10. (b). Local mean flow is from left to right in the 
horizontal direction. The direction of q in the plane is shown 
by arrows. The black contour corresponds to | q |= gth- 
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FIG. 4: Effective viscosity vs effective volume fraction <jf>e: red 
+ <p = 0.21, green x </> = 0.315. Daslied line Eilers fit HI with 
4>m = 0.58 and B = 1.5 [343 • Inset: relative increment of 
the volume fraction eis a function of Re. 



of a pair of spheres). 

We try to estimate the volume of the wake region by 
calculating the relative particle flux (relative momentum 
increments), defined as 



fi{r,r) = g(r,f)((5v}(r,f) 



(2) 



where (■) denote ensemble averaging and (5v is the rela- 
tive velocity of a pair of spheres. The relative particle 
flux in the shear plane is plotted in Figure [31 Clearly, 
the flux is largest in the region close to the surface of the 
sphere (i.e., grazing incidents) and at z > 3a (where the 
mean flow determines the flux). Most important, there 
exists a minimum, close to zero, behind a sphere when 
Re > 1. To estimate the volume occupied by the wake, 
we select a threshold value qth = 0.03 (black contour in 
figure [3]) and calculate the volume of the region where 
I q |< Qth- Denoting this volume, function of the par- 
ticle Reynolds number and volume fraction, Vd(Re, 0), 
the relative increase of excluded volume A(j)/(j) = Vd/Vg, 
where Vg = 47r(2a)''/3. The effective volume fraction, 
4>e = (fi + A(/)(Re), is plotted as a function of Re for two 
different values of 4> in the inset of Fig.|4l The increase of 
the volume fraction is significant, of the order of 10% for 
Re > 1. At fixed Re, the relative increase of 0e decreases 
marginally at larger (jj since collisions among particles are 
more frequent and deflect the particle trajectories reduc- 
ing the size of the wake region. The values of the effective 
viscosity, /i//io, in the range of cf) and Re considered can 
be collapsed to an universal curve using the effective vol- 
ume fraction (j)e((l>. R e), see FiglH where we also plot the 
Eilers fit (H]) jl4l4la |. valid for suspensions of vanishing 
inertia. We indeed find a good agreement between Eilers 
Fit and our data given the crude nature of estimate of 
the relative increase of volume fraction [i^ . 

We stress that the increase in effective volume frac- 
tion is essentially due to the formation of anistropic mi- 
crostructures, as already seen in Figs. [2] and [31 We 
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FIG. 5: (a) The pair-distribution function averaged over unit 
sphere, G{r) = (l/47r) / g{r,f)dQ, vs r/a - 2 for = 0.315 
and several values of Re: red -I- Re = 0.1, green x Re = 1., 
blue * Re = 5. and magenta □ Re = 10.. Inset: Second order 
structure function of longitudinal velocity differences of the 
spheres vs r/a — 2. (b) Probability distribution function of 
the number of clusters formed by A'^ spheres. 



present three further evidences to support this claim: 
(a) We plot in Fig. [5^ the pair-distribution function av- 
eraged over the solid angle G{r) = (l/47r) / g{r,f)dVL. 
Although clustering at small distance is clearly present 
(G(r) > 1 for small r), no significant change is ob- 
served with Re. (b) We report the second order structure 
function of the longitudinal particle velocity difference, 
S2{r) = (l/47r) (r)^)dr2, as an inset in figure [5^. 
5v\\ (r) = [vp — vq] • f , where Vp and vq are the veloci- 
ties of the P-th and Q-th particle separated by a distance 
r. Similar to G(r), 5*2 (r) does not show any significant 
change at small separation r when increasing Re. (c) We 
display in Fig. [SJd the probability distribution function 
of the number of clusters containing N spheres, Pc{N). 
Particles are considered to belong to the same cluster if 
their gap distance is less than 2% of a. We find that 
Pc{N) ^ N^'^, i.e., there exists finite probability to find 
large aggregates, similar to shear-thickening colloidal sus- 
pensions [j, However Pc{N) does not change as a 
function of Re. Hence, though hydroclusters are present, 
we do not observe a direct connection between formation 
or increase of clusters and inertial shear-thickening. 



4 



0.5 
0.4 
0.3 

r 

y 
0.2 

0.1 



0.5 



0.4 









■ w \ 






ft 't If • / 
/# - 

■ /' y</ 








Re=0.1 






Re=1.0 


---X--- 


/ (a) 


Re=5.0 


■■■«■■■■ 


Re=10. 


{3 



0.1 



0.2 ^/jj 0.3 



0.4 



0.5 



0.3 



0.2 
0.1 



0.006 





0.2 0.3 

z/H 



FIG. 6: Wall-normal profile of (a) the average local volume 
fraction, and (b) average particle velocity, V{z), for 4> ~ 0.315. 
The inset shows K = v^^s versus the wall distance z/h. Sym- 
bols are as in FigO 



The wall-normal profile of the local mean volume frac- 
tion (p{z), Fig. [5^, shows that particles tend to form layers 
due to the confinement from the wall, see also the mean 
particle velocity V{z) in \Ejp. Again layering does not 
show a monotonic behavior with Re. Interestingly, the 
single-point particle velocity fluctuations decrease with 
the inertia, as shown in the inset of Fig. The system 
appears more stable, with a more ordered structure and 
fewer particles jumping among the layers. 

The decrease of the fluctuation level and the increase 
of ordering is consistent with the idea of an increasing 
effective volume fraction at high shear-rates: the system 
tends to freeze as there is less available space for the 
particle motion. We argue that when the effective vol- 
ume fraction approaches the critical packing, the system 
would jam with vanishing fluctuations. Hence we sug- 
gest that the discontinuous shear-thickening observed at 
high concentrations, higher than those simulated here, 
yet below the geometrical maximum packing , can be 
interpreted as an increase of the effective volume fraction 
above </* < (f'm < 0e(Re). This behavior might ap- 
pear as heterogeneity in space with part of the system 
jammed at large shear rates 7 [lol |. 

To conclude, we show in this letter that inertial shear 



thickening in non-Brownian suspensions can be under- 
stood in terms of an increase in the effective volume 
fraction of the suspension. This, in turn, is due to for- 
mation of anistropic microstructures. We hope our work 
initiates an experimental search for such microstructures 
using e.g., magnetic resonance imaging. 
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